function energy=US_energy(dados,t_axis,y_axis)
tamanho=size(t_axis,2);
energy=zeros(1,tamanho);

t=(1:tamanho)/40e6;
inicio=t_axis(1)
fim=t_axis(tamanho)
iter=1;
for i=inicio:fim
    energy(iter)=trapz(dados(y_axis,i).^2);
    iter=iter+1;
end

figure
plot(energy)
grid on
xlabel('tempo (s)');ylabel('energia');
title('energia')